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ABSTRACT 

MARVIN  MARCUS:  Eigenvalue  localization; 
quadratic  and  bilinear  numerical  ranges; 
symmetry  classes  of  tensors 

Research  -eii^her  appeared  or  in  press  during  the  period 
October  1,  19§TT  "September  3.0 y _19^  covered  the  following 

topics:  (1)  generalizations  of  the  Cauchy-Schwar z  inequality 

to  determinants;  (2)  the  structure  of  the  set  of  complex 
numbers  (Axj  +  (Axr  ,yx  )  where  the  inner  products 

)/  {Ji  ,Jj  )  and  are  stipulated,  i  =  l,...,r; 

(3)  the  positive  definite  property  of  the  A-numerical  radius 
defined  by  Goldberg  and  Straus;  (4)  conditions  for  two 
decomposable  tensors  in  a  symmetry  class  to  be  equal  or  for 
a  decomposable  tensor  to  be  0;  (5)  the  geometry  of  the  set 
of  numbers  (Ax,y)  in  which  x  and  y  are  unit  vectors  with 
fixed  inner  product.  ___ 

This  report  summarizes  the  current  research  of  M.  Marcus 
sponsored  by  the  Air  Force  Office  of  Scientific  Research 
covering  the  general  topic  of  the  geometry  of  the  quadratic 
and  bilinear  numerical 'ranges ► 

1.  Generalizations  of  the  Cauchy-Schwarz  inequality  to 
determinants. 

In  answer  to  a  research  problem  posed  by  A.  Abian  [1] , 
we  were  led  to  consider  the  inequality  (1)  defined  below. 

Let  A,  B,  P,  Q  be  n-square  complex  matrices  and  for 
Ism  <  n,  let  U  and  V  be  n  x  m  complex  matrices.  We 


1 


2 


examine  the  relations  which  must  exist  between  A,  B,  P,  and 
Q  so  that 


det  (U*  AV)  det  (V*BU)  S  det  (U*PU) det (V *QV)  (1) 

holds  for  all  U  and  V.  If  m  =  1  and  A,  B,  P,  Q  are  taken  to 
be  the  n-square  identity  matrix,  then  (1)  reduces  to  the 
Cauchy-Schwar z  inequality 

2 

|(u,v)  I  S  { u , U )  ( V  , V )  , 

where  u  and  v  are  arbitrary  n  *  1  column  vectors. 

If  H  is  hermitian  positive  definite  (negative 
definite) ,  we  write  H  >  0  (H  <  0) .  If  H  and  K  are  definite 
of  the  same  sign  (opposite  signs) ,  we  write  HK  >  0  (HK  <  0) . 
The  principal  result  generalizing  the  Cauchy-Schwar z 
inequality  is  contained  in  the  following  theorem  [13]. 

Theorem.  Assume  n  S  3,  1  r_;  m  <  n,  and  A,  B,  P,  0  are 
nonsingular  complex  matrices.  Then  (1)  holds  for  all  n  x  m 
rectangular  matrices  U  *and  V  if  and  only  if 

(i)  P  =  aH,  Q  =  0K,  (a0  )m  =  e  =  ±1,  H  and  K  are 

definite  hermitian  and 

(ii)  A*  =  o)B,  uj  €  C  ,  uum  =  X€F,so  that  (1)  reads 

X  |  det(V*BU)  | 2  s  e  det (U*HU) det (V*KV) 


(2) 
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and  if 


are  the  eigenvalues  of  H 

1  n 


-1 


B*K-1B  and 


£  = 


a 


x  ; 

max  —  H 

Q  C  i=l 

m,n 


Ca(i) 


then: 


(iii)  if 
and  £  Si,  or  m 
or 

(iv)  if  e 
or  m  is  even  and 

(v)  if  c  = 


e  =1,  X  >  0,  then  either  m  is  odd,  HK  >  0 , 
is  even,  HK  can  have  either  sign,  and  £  -  1; 

=  1,  X  <  0,  then  either  m  is  odd  and  HK  >  0, 
HK  can  have  either  sign;  or 

-1,  X  >  0,  then  m  is  odd,  HK  <  0,  and  £  §  1; 


or 

(vi)  if  e 


-1,  X  <  0,  then  m  is  odd  and  HK  <  0. 


We  include  a  table  of  conditions  (iii)  -  (vi) : 


E 

X  . 

.  m 

H 

K 

£  . 

1 

-1- 

odd 

+ 

+ 

5  1 

1 

+ 

even 

^  1 

1 

-  . 

odd 

+ 

+ 

1 

- 

even 

-1 

+ 

odd 

+ 

5  1 

-1 

- 

odd 

+ 

It  should  be  noted  that  for  n  =  2,  m  =  1  the  inequality 
(2)  can  hold  with  both  H  and  K  indefinite.  Simply  take 


A 
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[i  0 

B  =  V  H  -  K  -  [o  -1 

2.  The  G-bilinear  range . 

The  numerical  range  of  a  linear  operator  A:V  ->  V, 
denoted  by  W(A),  is  the  image  of  the  surface  of  the  unit 
ball  under  the  quadratic  mapping 

x  +  (Ax,x) .  (3) 

That  is, 

W  (A)  =  {  (Ax , x)|  l|x|l  =  1)  . 

In  1918  and  1919  [2,3]  Toeplitz  and  Hausdorff  proved  the 
celebrated  result: 


W ( A )  in  a  convex  set  in  the  plane. 

Since  that  time  a  long  sequence  of  proofs,  generalizations 
and  refinements  of  this  result  have  appeared  [4] .  Perhaps 
the  two  most  significant  of  these  appeared  in  1937  and  1951. 
The  first  of  these  was  a  paper  by  von  Neumann  [5]  in  which 
the  values  of 


tr (XAYB) 


(4) 
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were  examined  for  A  and  B  fixed  linear  operators  on  V  with  X 
and  Y  varying  over  all  unitary  operators.  von  Neumann 
showed  that  the  totality  of  values  (4)  is  a  disk  centered  at 
the  origin  whose  radius  is 


n 


Z 

i=l 


(5) 


where  a.  g  •  •  •  s  a  and  3,  a  •••  §  3  are  the  singular 

1  n  1  n 

values  of  A  and  B  respectively.  The  second  result,  due  to 
Fan  [6] ,  dealt  with  values  of 


xi) 


(6) 


in  which  Xj,...,x  are  arbitrary  orthonormal  (o.n.)  vectors. 
Of  course  (6)  is  similar  in  form  to  (4) .  In  fact,  the  set 
of  values  (6)  is  precisely  the  same  as  the  set 

Wr(A)  =  {tr(PA)|  P  is  an  orthogonal  projection  of  rank  r}.  (7) 

» 

In  his  book  [7]  Halmos  calls  the  set  Wr (A)  the  rth  higher 
numerical  range.  This  set  is  also  convex  for  any  A,  a  fact 
first  proved  by  Berger  [8]  for  r  §  1.  More  recently, 
Goldberg  and  Straus  [9]  considered  the  C-numerical  range  of 
A: V  ->  V 


WC(A)  = 


{tr (CU^AU)  |  U  unitary} 


(8) 


6 


where  C  is  a  fixed  operator  on  V.  It  is  not  true  that  (8) 
is  a  convex  set  for  any  C,  however  if  C  is  hermitian  then 
WC(A)  is  convex  for  any  A,  a  result  originally  due  to  Westwick 
[10]  and  quite  recently  reproved  using  algebraic  methods  by 
Poon  [11]  ,  ' 

The  numerical  range  can  be  regarded  as  a  special  case 
of  the  q-bilinear  range 

W  ( A :  q )  =  {(Ax  ,y)  |  ]|xj|  =  |jy||  =  1,  (x,y)  =  q}.  (9) 

If  q  =  1,  W(A:q)  =  W (A) .  In  1961  Bush  and  Olkin  [12]  showed 
that  if  V  is  a  Euclidean  space  over  P,  q  >  0,  and  A  is 
symmetric  then  W (A:q)  is  the  real  interval 


[((1 

+  q)Xn 

-  (1 

-  q)A1  )/2  ,  (  (1  +  q)  Xx  - 

(1  - 

-  q) x  )/2] r 

n 

(10) 

where 

X1  * 

•  •  *  ^ 

are  the  eigenvalues 

pf 

A.  In  [9] 

the 

resul  t. 

(10)  i 

s  gene 

rali zed  In  various  ways. 

For  1 

s  r 

$  n  let  G  be  an  r  x 

3r 

complex  mal 

tr  ix 

partitioned  into  three  ,r  x  r  submatrices: 

G  =  [G1:G2:G3J  .  (11) 

A  set  of  2r  vectors  x, , . . .  ,x  ,  y,  » • • • »V  are  G-vectors  if 

1  r  1  J  r  - 

[  (xirXj)  1  =  Gj 


f 


(12) 
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t  (Y±  fYj  )  1  =  G  2' 


(13) 


and 


f  (x±  ,y  )  1  =  g3- 


(14) 


We  shall 

definite 

G-vectors 

y,*--- »yr  * 

function 


assume  henceforth  that  Gj  and  G2  are  positive 
hermitian  so  that  if  x  ,...,xr  and  y1,.../yr  are 
then  Xj,...rx  are  linearly  independent  as  are 
In  this  paper  we  investigate  the  values  of  the 


r 

cp(x;y)  =  2  (Axk  ,  yk) 

k=l 


(15) 


as  x,  ,...rx  and  v.  ,...,y  run  over  all  sets  of  G-vectors. 
1  r  *  1  r 

We  define  the  G-bilinear  range  to  be  the  set 


VI  ( A ;  G ) 


:  -  <  ?  y  ) 


.  ,v  G-vectors).  (16) 
"*  r 


Of'  course  the  set  W(A;.G)  may  be  empty  for  certain  choices  of 
G.  In  fact  we  shall  say  that  G  is  admi ssable  if  and  only  if 
there  exist  G-vectors.  Henceforth  we  shall  assume  that  G  is 
admissable  so  that  W(A;G)  is  not  empty.  In  Theorem  3  we 
will  obtain  some  necessary  and  sufficient  conditions  for 
admissability. 

The  principal  results  appearing  in  [  14 1  are 
Theorem  1.  If  1  §  r  s  n,  gig2  =  G2Gl'  an^ 
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G  =  Gly^2  G1^2 
3  1  2 

then  for  any  A,  W(A;G)  is  a  convex  set 

Theorem  2.  If  G0  f-  XI  and  tr  (Gn  ) 

—  J  n  -  3 

origin  i f  a'nd  only  if  A  =  0 . 

Theorem  2  is  closely  related  to  a 
Straus  [9] .  Note  that  if  r  <  n  then 
automatically  holds. 

The  conditions  for  admissability  are: 

Theorem  3.  The  r  x  3r  matrix  G  =  [G^iG^rG^]  is  admissable 
for  the  n-dimensional  unitary  space  V  i^f  and  only  i f 

a)  the  largest  eigenvalue  of 

G11g3G21(33*  08) 

does  not  exceed  1  and 

b)  (18)  has  1  as  an  eigenvalue  its  multiplicity  k 

s  a  t  i  s  f  j  e  s 


(17) 

in  the  plane. 

f-  0  then  VI  ( A ;  G )  j^s  the 

result  of  Goldberg  and 

the  conditions  G  f  XI 
3  n 


.  n  -  r  a  r  -  k. 

3.  Elementary  proofs  of  the  Goldberg-S  tr  a us  theorem . 

For  a  pair  of  n-square,  complex  matrices  A  and  Bf  M. 
Goldberg  and  E.G.  Straus  [91  have  defined  the  A-numerical 
radius  of  B  as  the  quantity 


rA  (B)  =  max  |  tr  (AU*BU)|  , 
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where  the  maximum  is  taken  over  all  unitary  matrices  U.  In 
a  later  paper  [ 15 1  they  showed  that  rA  is  a  generalized 
matrix  norm  if  and  only  if  the  conditions 

A  is  nonscalar  and  tr  A  4  0  (19) 

hold.  A  generalized  matrix  norm  is  a  scalar  valued, 
homogeneous,  subadditive,  positive  definite  matrix  function. 
The  homogeneity  and  subadditivity  properties  are  easily 
established,  as  is  the  necessity  of  the  condition  (19)  from 
the  positive  definiteness  of  rA  .  Thus  the  majority  of 
section  2  in  the  previously  cited  paper  is  concerned  with 
showing  that  condition  (19)  implies  the  positive 
definiteness  of  the  A-numerical  radius. 

While  examining  the  proof  of  this  last  statement  in  a 
seminar  held  during  the  summer  of  1980  in  which  M.  Goldberg 
participated,  three  separate  proofs  of  this  result  were 
obtained  by  the  present  authors. 


4.  Symmetry  classes  of  tensors . 

In  a  paper  by  M.  Marcus  and  J.  Chollet  [17  ],  necessary 
and  sufficient  conditions  were  stated  for  two  decomposable 
symmetrized  tensors  to  be  equal.  There,  the  linear 
independence  of  the  vectors  forming  such  tensors  was 
assumed.  In  [  24]  this  assumption  is  dropped  and  more 


general  requirements  for  equality  are  obtained.  Conditions 
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for  a  decomposable  symmetrized  tensor  to  be  zero  are  also 
given. 

A  brief  history  of  the  above  problem  is  as  follows. 
Before  the  paper  of  Dias  da  Silva  [18  ]  ,  most  results 
concerned  symmetry  classes  whose  associated  characters  were 
of  degree  "one  (see  theorem  4.5,  p.136,  of  [19]).  R.  Merris 
[20]  did  consider  more  general  symmetrizers  and  higher 
degree  characters  but  did  not  give  necessary  and  sufficient 
conditions  for  equality.  This  was  done  by  Dias  da  Silva  in 
[18]  and  then  with  de  Oliveira  in  [21].  A  different 
approach  was  initiated  in  [17]  and  is  continued  here. 
Related  results  have  appeared  in  [22]  and  [ 23 ] - 


5.  The  q- bilinear  range 

Numerical  experiments  are  currently  under  way  on  the 
q-bilinear  range 

W(A:q)  =  { ( Ax ,  y )  j  (x,x)  =  (y,y)  =  1,  (>:,”)  =  cn. 

It  is  conjectured  that  the  set  W(A:q)  is  convex  for 
0  <  q  <  1.  For  n  =  2,  the  following  graphic  outputs  were 
obtained  using  an  APPLE  II-plus  computer  and  the  GRAPHTRIX 
plotting  software  by  DATA  TRANSFORMS,  INC. 


1.  A. 
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I,  Marcel  Dekker,  New  York,  1973. 

20.  R.  Merfis,  Equality  of  Decomposable  Symmetrized 
Tensors,  Canad.  J.  Math.  27(1975),  1022-1024. 

21.  G.N.  de  Oliveira  and  J.A.  Dias  da  Silva,  Conditions  for 
Equality  of  Decomposable  Symmetrized  Tensors  II,  Linear 
Algebra  and  Appl.  28(1979),  161-176. 

22.  M.  Marcus,  Decomposable  Symmetrized  Tensors  and  an 
Extended  LR  Decomposition  Theorem,  Linear  and  Multilinear 
Algebra,  6(1978),  327-330. 

23.  M.  Marcus  and  J.  Chollet,  Linear  Groups  defined  by 
Decomposable  Tensor  Eaualities,  Linear  and  Multilinear 
Algebra  8(1979/80),  20*7-212. 

24.  M.  Marcus  and  J.  Chollet,  On  equality  of  decomposable 
symmetrized  tensors,  in  press,  Linear  and  Multilinear 
Algebra . 

25.  M.  Marcus  and  J.  Chollet,  The  index  of  a  symmetry  class 
of  tensors,  in  press.  Linear  and  Multilinear  Algebra. 
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Marvin  Marcus 

Publications:  October  1,  1980  to  date 

1.  M.  Marcus  and  K.  Moore,  A  determinant  formulation  of  the 
Cauchy-Schwar z  inequality,  Linear  Algebra  Appl.  36(1981), 
111-127. 

2.  M.  Marcus  and  J.  Chollet,  On  the  equality  of 
decomposable  symmetrized  tensors,  to  appear.  Linear  and 
Multilinear  Algebra. 

3.  M.  Marcus  and  J.  Chollet,  The  index  of  a  symmetry  class 
of  tensors,  to  appear,  Linear  and  Multilinear  Algebra. 

4.  M.  Marcus  and  M.  Sandy,  Three  elementary  proofs  of  the 
Gold^rg-Straus  theorem  on  numerical  radii,  to  appear. 
Linear  and  Multilinear  Algebra,  11(3),  1982. 

5.  M.  Marcus,  K.  Kidman  and  M.  Sandy,  Bilinear  ranges,  in 
preparation. 

6.  Solution  to  Problem  #1122  (with  M.  Newman),  Mathematics 
Magazine,  to  appear. 


Abstract 


M.  Newman:  THE  APPLICATIONS  OF  NUMBER  THEORY  TO  COMPUTATION 

Research  completed  in  the  last  period  includes  the 
following : 

(1)  The  preparation  of  programs  using  modular  arithmetic  for 
serial  machines  to  find  all  rational  solutions  of  a  system  of 
linear  equations.  (2)  The  preparation  of  programs  using 
multiprecision  integral  arithmetic  for  serial  machines  to  find 
the  Hermite  normal  form  and  the  Smith  normal  form  of  an 
integral  matrix.  (3)  The  preparation  of  programs  using 
modular  arithmetic  for  the  computation  of  the  elementary 
divisors  of  an  integral  matrix. 
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M.  NEWMAN:  The  Applications  of  Number  Theory  to  Computation 

The  purpose  of  this  report  is  to  summarize  my  Air  Force 
sponsored  research  on  the  applications  of  number  theory  (in 
particular  modular  arithmetic  and  multiprecision  integral 
arithmetic)  to  problems  of  computation,  using  available  serial 
computers  (ILLIAC  IV  unfortunately  no  longer  exists) . 

The  principal  objective  of  my  research  is  to  apply 
number-theoretic  ideas  to  problems  of  numerical  analysis  and 
computation.  In  particular,  applications  of  modular  and 
multiprecision  arithmetic  are  being  made  to  the  following 
topics : 


(1)  The  exact  solution  of  an  integral  system  of  linear 
equations,  and  the  exact  computation  of  the  determinant 
of  the  system; 


(2) 

The 

determination 

of 

the 

exact 

inverse  of 

an 

integral  matrix  using  minimal 

storage ; 

(3) 

The 

determination 

of 

the 

rank  and 

a  basis  for 

the 

n  u  1 1 

space 

of  an  in  teg  r  co¬ 

ma  t  r i x ; 

(4) 

The 

determination 

of 

all 

rational 

solutions  of 

an 

integral  system  of  -linear  equations; 

(5)  The  computation  of  the  Hermite  normal  form  of  an 

integral  matrix; 

(6)  The  computation  of  the  Smith  normal  form  of  an 

integral  matrix; 

(7)  The  computation  of  the  elementary  divisors  of  an 

integral  matrix. 

Programs  to  perform  (1) ,  (2) ,  (3) ,  and  (4)  using  modular 
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arithmetic  have  been  prepared  and  checked.  The  programs  are 
written  in  a  highly  portable  version  of  FORTRAN  and  were 
implemented  without  difficulty  on  a  number  of  different 
computers,  including  a  microcomputer. 

Programs  to  perform  (5) ,  (6) ,  and  (7)  have  been  prepared 
in  two  versions:  the  first  set  of  programs  uses  the  standard 
integer  arithmetic  provided  by  the  host  computer,  and  are 
limited  by  the  usual  problems  of  integer  overflow.  The  second 
set  uses  multipc ecision  integer  arithmetic  subroutines  written 
by  the  proposer,  and  are  guaranteed  to  produce  the  correct 
result.  Und.=  investigation  is  the  problem  of  whether  or  not 
residue  arithmetic  may  be  employed  in  the  computation  of  these 
important  normal  forms. 

Two  Master^s  Theses  on  this  topic  were  written  under  my 
supervision  during  this  period,  by  G.  Rasmussen  and  R.  Demb. 

A  report  on  the  work  performed  is  under  preparation  and 
will  appear  in  Math.  Comp,  at  some  future  time. 


J 


M.  Newman 


PUBLICATIONS  October  1,  1980  to  date 

1.  Matrices  of  finite  period  and  some  special  linear 
equations.  Linear  andd  Multilinear  Algebra  8  (1980) ,  189-195. 

2.  On  a  problem  suggested  by  Olga  Taussky-Todd ,  Illinois  J. 
Math.  24  (1980),  156-158. 

3.  A  surprising  determinantal  inequality  for  real  matrices 
(with  C.R.  Johnson),  Math.  Ann.  247  (1980),  179-186. 

4.  A  note  on  cospectral  graphs  (with  C.R.  Johnson),  J.  Comb. 
Theory  B  28  (1980),  96-103. 

5.  Positive  definite  matrices  and  Catalan  numbers  (with  F.T. 
Leighton),  Proc.  Araer.  Math.  Soc.  79  (1980),  177-181. 

6.  Gersgorin  revisited.  Linear  Algebra  and  its  Applications 
30  (1980) ,  247-249. 

7.  Determinants  of  abelian  group  matrices  (with  M.  Mahoney), 
Linear  and  Multilinear  Algebra  9,  (1980)  ,  121-132. 

8.  Determinants  of  circulants  of  prime  power  order,  Linear 
and  Multilinear  Algebra  9  (1980),  187-191. 

9.  Simple  groups  of  square  order  and  an  interesting  sequence 
of  primes  (with  D.  Shanks  and  H.C.  Williams) ,  Acta  Arith.  38 
(1980) ,  129-140. 

10.  A  radical  diophantine  equation,  to  appear  in  J.  Number 
Theory. 

11.  Cyclotomic  units  and  Hilbert^s  Satz  90,  to  appear  in  Acta 

Arith . 

12.  Pseudo-similarity  and  partial  unit  regularity  (with  F.J. 
Hall,  R.E.  Hartwig,  and  I.J.  Katz),  to  appear  in  Czech.  J. 
Math. 

13.  Similarity  over  SL(n,F) ,  to  appear  in  Linear  and 
Multilinear  Algebra. 

14.  On  a  problem  of  H.J.  Ryser,  to  appear  in  Linear  and 
Multilinear  Algebra. 

15.  Lyapunov  revisited,  to  appear  in  Linear  Algebra  and  its 
Applications. 


16.  A  result  about  determinantal  divisors,  to  appear  in  Linear 
and  Multilinear  Algebra. 
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Abstract 

R.C .  Thompson:  MATRIX  EQUALITIES  AND  INEQUALITIES 
An  ongoing  program  studying  singular  values  and  invariant 
factors  has  been  continued.  Recent  work  includes  a  partially 
successful  attempt  to  analyze  how  the  invariant  factors  of 
integral  matrices  behave  when  matrices  add.  Nine  papers  have 
been  submitted  to  journals  {four  so  far  accepted) ,  with  two 
more  nearly  ready  for  submission,  and  several  others  being 
prepared.  Invitations  have  been  accepted  to  write  the 

Encyclopedia  of  Mathematics  volume  on  Matrix  Inequalities,  and 
to  deliver  the  keynote  lecture  on  "Core  Linear  Algebra"  at  the 
April  1982  SIAM  meeting  on  Applied  Linear  Algebra. 
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R.C.  THOMPSON:  Matrix  Equalities  and  Inequalities 
The  proposed  research  was  organized  into  four 
compartments.  I  discuss  these  in  turn. 

(1)  _ The  matrix  exponential.  I  have  conjectured  that  if  H,K 

are  Hermitian  matrices,  unitary  matrices  U,V  will  always  exist 
so  that 

iH  iK  i  (UHU*  +  VKV*)  i  =  J~1  ,  e  =  2.718... 

g  g  =  6 

*  *  denotes  conjugate  transpose 

where  U,V  depend  on  H,K.  A  renewed  attempt  was  made  to  prove 
this  conjecture,  but  success  was  minimal.  Numerous  special 
cases  can  be  settled,  and  the  key  concepts  that  are  needed 
from  Lie  theory  have  been  identified,  but  the  unravelling  of 
the  details  remains  beyond  reach.  The  author's  belief  is  that 
several  years  of  sustained  effort  will  be  needed  to  settle 
this  conjecture. 

Several  manuscripts  bearing  on  this  question  have  been 
prepared,  but  none  have  been  submitted  to  journals. 

(2) ,  (3)  Matrix  valued  inequalities.  If  A  is  a  matrix,  let 

,  ,  *  ); 

I  A |  --  (AA  )~  be  its  matrix  valued  absolute  value.  T  proved 

several  years  ago  that  if  A, B  are  matrices,  then 

|a  +  b|  §  u  |  a  |  u*  +  v  |  b  |  v* 

where  unitaries  U,V  depend  on  A,B.  Here  g  means  the 
difference  of  the  two  sides  is  semidef inite.  It  has  now  been 
shown  that  this  inequality  is  valid  if  the  matrices  have 
quaternion  entries.  Whether  this  is  so  was  specifically  asked 
in  the  proposal,  and  now  has  been  settled  affirmatively. 
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Several  years  ago  I  conjectured  that  if  A  is  a  matrix 
with  U  its  unitary  factor  from  the  polar  decomposition,  and  if 
V  is  any  other  unitary,  we  always  will  have 

|  A  -  U  I  §  %  WL  |  A  -  V  j  WL  *  +  >5  W2  j  A  -  V  |  W  * 

for  appropriate  unitaries  W  ,  W  .  That  is:  U  is  nearer  A  than 
any  other  unitary  V  is.  This  conjecture  is  true  for  2x2 
matrices,  but  cannot  so  far  be  settled  for  the  3x3  case,  in 
spite  of  copious  quantities  of  computer  work  that  repeatedly 
verifies  its  truth  in  numerical  cases. 

Considerable  (theoretical)  effort  has  been  expended  on 
this  conjecture,  with  not  much  success.  Further  work  is 
planned.  This  is  one  of  the  items,  not  explicitly  listed  in 
the  grant  proposal,  but  caught  under  the  general  theme  of 
"study  of  matrix  valued  inequalities". 

(4) _ Harmonic  analysis.  A  matrix  question  coming  from 

harmonic  analysis,  previously  solved  by  this  author,  was 
written  into  a  paper  and  submitted  to  a  journal  (#4  below) . 

Other  work  not  mentioned  in  the  grant  proposal. 

(i)  How  do  the  invariant  factors  of  integral  matrices 
behave  when  the  matrices  add?  A  five  month  investigation 
resulted  in  a  complete  solution  for  2x2  and  3x3  matrices, 
and  a  conjectured  solution  for  n  x  n  matrices. 

(ii)  An  invitation  has  been  accepted  to  write  the  Matrix 
Inequalities  volume  for  the  F.ncyclopedia  of  Mathematics.  The 
Table  of  Contents  for  this  volume  is  attached  to  this  report. 
It  will  take  several  years'  work  to  bring  this  volume  to 


22 


completion.  Work  on  it  was  begun  under  the  just  expired 
grant. 

(iii)  An  invitation  has  been  accepted  to  deliver  one  of 
the  principal  lectures  at  the  SIAM  Applied  Linear  Algebra 
meeting,  April,  1982.  This  meeting  has  six  principal 
speakers,  five  on  various  aspects  of  applied  linear  algebra, 
plus  (to  keep  the  meeting  honest)  one  on  core  linear  algebra. 
I  am  the  invited  speaker  on  core  linear  algebra,  and  work  on 
it  was  begun  under  the  just  expired  grant. 


R.C.  Thompson 

PUBLICATIONS  October  1,  1980  to  date 

1.  P-adic  matrix  valued  inequalities,  submitted. 

2.  The  Jacobi-Gundelf inger-Frobenius-Iohvidov  rule  and  the 
Hasse  symbol,  to  appear  in  Letters  in  Linear  Algebra. 

3.  The  Smith  form,  the  inversion  rule  for  2x2  matrices,  and 
the  uniqueness  of  the  invariant  factors  for  finitely  generated 
modules,  to  appear  in  Letters  in  Linear  Algebra. 

4.  Doubly  stochastic,  unitary,  unimodular,  and  complex 
orthogonal  power  embeddings,  submitted. 

5.  A  converse  to  the  Ostrowski-Taussky  determinantal 
inequality,  to  appear  in  Linear  Algebra  and  Applications. 

6.  Cyclic  relations  and  the  Goldberg  coefficients  in  the 
Campbell-Baker-Hausdor f f  formula,  submitted. 

7.  An  inequality  for  invariant  factors,  submitted. 

8.  The  matrix  valued  triangle  inequality  -  quaternion 
version,  to  appear  in  Linear  and  Multilinear  Algebra. 

9.  The  true  growth  rate  and  the  inflation  balancing 
principle,  submitted. 

The  following  two  manuscripts  are  essentially  complete  and 
about  to  be  sent  to  journals. 

10.  Left  multiples  and  right  divisors  of  integral  matrices. 

11.  Sums  of  integral  matrices. 

Some  of  these  papers  represent  work  from  ear"  ■  t  r\?'  • 

grants,  some  from  this  grant  alone. 

Seven  manuscripts  are  in  various  stages  of  preparation, 
and  will  be  listed  in  future  grants  as  they  eventually  get  to 
journals . 
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A  Beginner's  Survey  of  Matrix  Inequalities 
by  R.C.  Thompson 

Table  of  Contents 

1.  Dominance  and  Lie  theory 

2.  Inequalities  for  matrices  with  continuous  entries 

3.  Inequalities  for  matrices  with  discrete  entries 

4.  Matrix  valued  inequalities 

5.  The  matrix  exponential 

6.  Decomposition  theorems 

7.  The  numerical  range 

8.  Nonnegative  matrices 

9.  Scalar  norms 

10.  Other  topics 


Abstract 

H.  Mine:  VAN  DER  WAERDEN  PERMANENT  CONJECTURE. 

THEORY  OF  PERMANENTS  AND  ITS  APPLICATIONS. 
INEQUALITIES. 

Research  during  the  period  included: 

(1)  Study  of  Egorycev's  proof  of  the  van  der  V7aerden 
conjecture.  The  proof  was  reformulated  and  somewhat  extended. 
Alternative  proofs,  different  in  parts  from  Egory^ev's  proof, 
were  devised. 

(2)  The  developments  in  the  theory  of  permanents  in  the  last 
four  years  were  surveyed  and  discussed.  Present  status  of  30 
conjectures  and  13  problems  was  discussed.  A  comprehensive 
bibliography  was  compiled. 

(3)  New  rearrangement  inequalities  for  products  and  sums  of 


nonnegative  numbers  were  obtained. 
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H.  MINC:  (I)  van  der  Waerden  Permanent  Conjecture 

The  most  famous  and  one  of  the  most  important  problems  in 
the  theory  of  permanents  is  the  so-called  van  der  Waerden 
conjecture:  if  A  is  a  doubly  stochastic  n  x  n  matrix,  then 

per (A)  §  n!/n  , 

where  equality  can  hold  if  and  only  if  all  the  entries  of  A 
are  1/n.  The  problem  actually  originated  in  a  theorem  of 
Konig  (1916)  who  proved  that  the  permanent  of  a  doubly 
stochastic  matrix  is  always  positive.  For  over  half  a  century 
the  conjecture  resisted  all  attempts  of  solution,  until  late 
in  1980  when  Egorycev  proved  it  completely.  One  of  the 
effects  of  this  historical  proof  was  that  much  of  the  research 
in  progress  on  permanents  of  doubly  stochastic  matrices, 
including  applications  to  the  dimer  problem,  enumerations  of 
various  combinatorial  problems,  etc.,  had  to  be  reevaluated. 

In  [5]  a  version  of  Egorycev's  proof  is  given  in  which 
some  of  the  auxiliary  results  are  somewhat  generalized.  In 

#  \S  . 

[4]  it  is  shown  that  the  concluding  argument  m  Egorycev  s 
proof  can  be  replaced  by  either  of  two  alternative  proofs. 

(II)  Theory  of  Permanents  and  Its  Applications 
The  last  four  years  witnessed  an  increased  research 
activity  in  the  theory  of  permanents  and  its  applications.  At 
least  75  new  research  papers  on  the  subject  and  similar 
publications  were  written  or  were  actually  published  during 
this  four-year  period;  this  represents  nearly  20%  of  the  total 
literature  on  permanents  in  the  last  170  years. 

In  (5]  the  developments  in  the  area  are  surveyed.  The 


paper  (56  typed  pages)  includes  a  detailed  discussion  of 


Egorycev's  proof  of  the  van  der  Waerden  conjecture  (see 
above),  a  section  on  bounds  of  permanents  of  (0 , 1) -matrices , 
and  on  their  applications  to  the  dimer  problem,  to  the  problem 
of  enumeration  of  Latin  rectangles,  and  other  topics.  The 
paper  contains  a  report  on  the  current  status  of  each  of  the 
20  conjectures  and  10  problems,  which  were  listed  as  unsolved 
in  Minces  "Permanents"  (1978) .  Additional  10  conjectures  and 
3  problems  are  proposed.  An  extensive,  virtually  complete, 
bibliography  is  appended,  each  item  with  a  short  abstract. 
Finally,  30  additional  references  (with  abstract)  are  listed 
for  the  period  preceding  1978;  these  failed  to  be  included  in 
the  bibliography  of  "Permanents". 

(Ill)  Inequalities 

Let  (a,r  a2j ami)  be  non-negative  m(-tuples,  j  =  1. n,  where  m,  g . . .  & 

=  I ,  and  ,  m,  =  M.  Let 

a'I1g«5IS...sQ^ilga,'JSaS2S...Sa^,S...Sa^Sa5nS...sQ^ 

be  the  M  numbers  ait,  i=l . m(,  /=!, — ,  n,  arranged  in  non-decreasing 

order,  and  " 

a*  a  a?,S  . . .  a* f£  or* .  £  o?,  fc . . .  *  a*  .,  g  . . .  £  a *, >  ft*.  * . . .  5:  a^n 

be  the  same  numbers  arranged  in  non-increasing  order. 

In  [1]  we  prove  inequalities  of  the  following  types.  If 

as  1,  i  =  1,...,  m. ,  j  =  1,...,  n,  then 
i-J  j 


I  n  n 

i*»  i  »-  i  j-i  i 


...  rt  (1  +  •  •  •  am„)  s  n  (1  -  •  •  •  <„>. 


If  i  1  f  •  •  •  in  j  f  D  *"  1  f  • » •  i  n  f  thsn 


i-1 i-l 


j-l .-I 


and 

fl  (1  +  atia2i . . .  a„j)  S  fl  (1  +  a*  a*, . . .  a*J. 

i-t  i-i 

In  either  case 

n  n  m, 

n  i  ^=n  z  «*• 

i-i *-i  i-i  i-i 

i 

Some  of  these  inequalities  were  used  to  establish  bounds  for 
permaments  of  nonnegative  matrices. 


H.  Mine 


PUBLICATIONS  October  1,  1980  to  date 

1.  Rearrangement  inequalities,  Proc.  Roy.  Soc.  Edinburgh  86A 
(1980),  339-345. 

2.  Inverse  elementary  divisor  problem  for  doubly  stochastic 
matrices.  Linear  and  Multilinear  Algebra  (to  appear). 

3.  Inverse  elementary  divisor  problem  for  nonnegative 
matrices,  Proc.  Amer.  Math.  Soc.  (to  appear). 

y/ 

4.  A  note  on  Egorycev^s  proof  of  the  van  der  Waerden 
conjecture.  Linear  and  Multilinear  Algebra  (to  appear) . 


5.  The  theory  of  permanents  1878-1981,  (submitted). 
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E.C.  JOHNSEN:  Elementary  Factorizations  of 
Stochastic  and  Doubly  Stochastic  Matrices 
E.C.  Johnsen  has  been  investigating  the  questions  of  (1) 
when  stochastic  matrices  factor  into  elementary  stochastic 
matrices  and  (2)  when  doubly  stochastic  matrices  factor  into 
elementary  doubly  stochastic  matrices.  Results  obtained 
earlier  on  problem  (1)  are  being  reworked  and  rewritten.* 

These  results  give  rather  general  sufficient  conditions  for 
the  elementary  factorization  of  stochastic  matrices.  Problem 
(2)  appears  to  be  considerably  more  difficult  than  problem 
(1) .  Partial  affirmative  results  on  problem  (2)  have  been 
obtained  by  the  author  for  the  case  of  orthostochastic 

I 

matrices  of  order  3.  It  is  known  that  not  every  doubly 
stochastic  matrix  of  order  n  §  3  and  that  not  every 
orthostochastic  matrix  of  order  n  s  4  (M.  Marcus)  can  be 

factored  into  elementary  doubly  stochastic  matrices.  Thus  the 
general  factorization  existence  question  is  still  open  for 
orthostochastic  matrices  of  order  3. 

These  problems  are  of  applied  interest  because  of  the 
direct  relation  between  elementary  factorizations  of 
stochastic  and  doubly  stochastic  matrices  and  sequential 
allocation  and  reallocation  .programs  for  multiple  commodities 
at  fixed  sources  and  destinations  (blending  problems) . 

I 

i 

*  Johnsen,  E.C.,  "Stochastic  Matrices  II.  Some  Factorizations 
into  Elementary  Stochastic  Matrices" 
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